Strain and field modulation in bilayer graphene band structure 
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Using an external electric field, one can modulate the bandgap of Bernal stacked bilayer graphene 
by breaking A — B symmetry. We analyze strain effects on the bilayer graphene using the extended 
Hiickel theory and find that reduced inter layer distance results in higher bandgap modulation, as 
expected. Furthermore, above about 2.5 A interlayer distance, the bandgap is direct, follows a 
convex relation to electric field and saturates to a value determined by the interlayer distance. 
However, below about 2.5 A, the bandgap is indirect, the trend becomes concave and a threshold 
electric field is observed, which also depends on the stacking distance. 

PACS numbers: 73.22.-f, 73.20.-r, 72.80.Rj 



Graphene is a two dimensional membrane of hexago- 
nally arranged carbon atoms. One obtains zero bandgap 
and linear dispersion around the Dirac point for a mono- 
layer of graphene p], [3, @, 0] • When two graphene lay- 
ers are stacked in Bernal (A — B) configuration, these 
monolayer features are lost and the dispersion becomes 
quadratic 0|. The bilayer configuration is very in- 
teresting for different applications, since the bandgap 
can be modulated from zero to few tens of an elec- 
tron volt by using an external out-of-plane electric field 
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This electric field breaks the A — B symmetry in the 
stacked bilayers and hence lifts the degeneracy, resulting 
in a bandgap opening. The bandgap depends on the cou- 
pling between the two graphene layers [H, @]. Therefore, 
one should be able to tune this bandgap by straining the 
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FIG. 2: Electronic structure of the valence band of the bilayer 
graphene. With decreasing stacking distance (c), new features 
appear in the electronic structure due to increased coupling 
between the two layers. The color bar is in eV. 
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FIG. 3: Electronic structure of the conduction band of the 
bilayer graphene with various stacking distances. 



FIG. 1: (color online) Electronic structure of bilayer 
graphene. E(k) diagrams are shown around the Dirac point 
with and without electric fields. For equilibrium stacking dis- 
tance c = 3.35 A, the dispersion is quadratic. With an ex- 
ternal electric field, a direct bandgap appears. For c = 3 A, 
bandgap modulation is larger due to increased hopping be- 
tween the two layers and remains direct. For c = 2 A, the 
bilayer becomes metallic with E = and has an indirect 
bandgap. 



bilayer. This is analogous to piezoresistivity, however the 
effect is enhanced exponentially due to bandgap modifi- 
cation. This could be very useful in sensors and other 
applications involving strain. Additionally, we find that 
by straining bilayer graphene, conduction and valence 
band dispersions also change. 

For electronic structure calculations, we use a semi- 
empirical extended Hiickel theory (EHT) with non- 
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FIG. 4: (color online) Bandgap modulation of strained bilayer 
with an external electric field. The dashed line is for equilib- 
rium c = 3.35 A, for which bandgap saturates at about 0.13 
eV. A convex relation is observed as a function of electric 
field above c~ 2. 5 A. By decreasing c, bandgap modulation 
increases due to increased interlayer hopping. 



FIG. 5: (color online) Bandgap modulation of strained bilayer 
with an external electric field. Below about 2.5 A stacking dis- 
tance, a concave relation is observed as a function of electric 
field, where bandgap opening occurs only above a threshold 
electric field (Et). Above Et, the bandgap linearly increases 
with external electric field. The inset shows variation of Et 
with c. 



orthogonal basis set. The detailed model has been de- 
scribed in Ref. The EHT parameters are transfer- 
able for different systems and have been benchmarked 
with the graphene band structure using generalized gra- 
dient approximation in the density function theory. This 
method is computationally very efficient and the non- 
orthogonality captures the bonding chemistry of vari- 
ous systems very well. It has been applied to armchair 
graphene nanoribbons [l2j], large Si systems for elec- 
tronic structure and quantum transport 13l [l4j and sin- 
gle molecule incoherent transport [15|, LL6[ utilizing mod- 
est computational resources. 

The C-C atomic distance is taken as 1.44 A within 
the graphene plane. Out-of-plane equilibrium stacking 
distance (c) is 3.35 A as shown in Fig. 1. The gray 
(light shaded) atoms belong to lower layer and the blue 
(dark shaded) atoms represent upper layer. The unit cell 
consists of four atoms - two atoms in each layer shown 
by red (black) dotted line in Fig. 1. The two gray (light 
shaded) atoms in the lower layer are referred to as A and 
B, respectively. The two blue (dark shaded) atoms in 
the upper layer are referred to as A and B, respectively. 
The electric field is applied in the out-of-plane direction 
as shown in Fig. 1. Atomic visualization is done using 
GaussView [17j. The substrate effects are ignored and 
in that sense, the graphene bilayer is assumed to be in 
vacuum. 

The electronic structure calculations for a bilayer 
graphene with equilibrium stacking distance c = 3.35 A 
is shown in Fig. 1. Without any electric field, the dis- 
persion is quadratic with a zero bandgap at the Dirac 
(K) point. With an electric field of 1 V/nm, a direct 
bandgap opens up due to A — B symmetry breaking and 
the conduction/ valence band minimum/maximum shifts 



away from the Dirac point. By decreasing c to 3 A, 
without an electric field, the bandgap is still zero with 
quadratic dispersion. However, the dispersion changes 
with an increase in the effective mass and the high lying 
bands move farther away from the Dirac point. With E 
= lV/nm, again a direct bandgap is observed with the 
same features as the ones for c = 3.35 A. Apart from 
this, the bandgap modulation is higher for c = 3 A, al- 
though the Laplace's potential is smaller for the same 
electric field due to the reduced spacing. This is due to 
the increased wave function overlap between the two lay- 
ers, which varies exponentially with distance as compared 
to the Laplace's potential which is linearly dependent on 
the distance. Further reducing the stacking distance re- 
sults in a metallic state as shown in Fig. 1 for c = 2 A. 
By applying 1 V/nm electric field, an indirect bandgap 
is created. 

Below about 2.5 A, the distance between atoms in the 
bottom layer and the next-nearest neighbors in the top 
layer becomes about 3.35 A, which is the equilibrium 
stacking distance for bilayer graphene. Therefore, the 
hopping integral between A — B sites becomes compa- 
rable to the equilibrium one of the A — B sites. This 
results in the reported electronic structure modifications 
in Fig. 1 for c = 2 A. This electronic structure change is 
also evident in the color plots of the valence band states 
over the two dimensional Brillouin zone in Fig. 2. With 
decreasing stacking distance from c = 3.35Atoc = 2.5 
A in the absence of an electric field, the features become 
broad indicating an increase in effective mass. Further 
decreasing c to 2 A results in directional change of the 
features, which signifies that the additional bonding hap- 
pens in different directions than that for c = 3.35 A. 
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These features are also present in the two dimensional 
Brillouin zone of the conduction band states in Fig. 3. 
However, it is more evident in the valence band states due 
to the constructive wave function overlap for the bond- 
ing nature in contrast to the anti-bonding nature of the 
conduction band states. 

Moreover, we find that the smaller the stacking dis- 
tance, the higher the bandgap modulation by electric 
fields as shown in Fig. 4. This is a desirable trait for prac- 
tical applications. The bandgap increases with increasing 
electric field and then saturates. For the equilibrium c 
= 3.35 A, the bandgap saturates at about 0.13 eV. For 
the local density approximation the bandgap satu- 
rates at about 0.25 eV. Although the results are different 
for the two techniques, the general trend after straining 
the bilayer should remain the same in the two methods. 
Above about 2.5 A, the trend of field modulation is con- 
vex and in the absence of electric field, a zero bandgap is 
observed as shown in Fig. 4. However, below about 2.5 
A, the systems becomes metallic without electric fields 
and the trend of field modulation becomes concave, for 
which a threshold behavior is observed as shown in Fig. 
5. The bandgap opening occurs only above a certain 
threshold electric field. This threshold value monotoni- 
cally increases with the decreasing stacking distance as 
shown in the inset of Fig. 5. As shown in Fig. 1, for 



c = 2 A, the valence and conduction bands are overlap- 
ping, giving rise to a metallic state. Therefore, a certain 
electric field needs to be applied to reduce this overlap 
to zero and then additional electric field opens up the 
bandgap. A bandgap of about 0.6 eV can be opened as 
shown in Figs. 3 and 4 using a 4 V/nm electric field. 
Calculations for such a high field are reported to show 
the bandgap saturation and the threshold trends for the 
strained bilayer. However, these high fields may not be 
feasible in devices due to the dielectric reliability concerns 
and physical constraints. High-K dielectrics may be used 
as an alternative to enhance field inside the graphene bi- 
layer, while keeping the dielectric within the breakdown 
regime. 

We have analyzed the strain effects on the bilayer 
graphene using EHT and the corresponding bandgap de- 
velopment due to an external electric field. We find that 
above 2.5 A stacking distance, the bandgap is direct and 
the bandgap follows a convex relation to the electric field. 
However, below 2.5 A stacking distance, the bandgap be- 
comes indirect and it follows a concave relation. 
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